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Abstract. In this note we will study the p-divisibility of multiple harmonic sums. In particular 
we provide some generalizations of the classical Wolstenholme's Theorem to both homogeneous and 
non-homogeneous sums. We make a few conjectures at the end of the paper and provide some very 
convincing evidence. 
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1 Introduction 

The Euler-Zagier multiple zeta functions of length I are nested generalizations of Riemann zeta 
function. They are defined as 

c(s) = c(si,...,so= E k-^'---kr' (1) 

0<fci<-<fei 

for complex variables si, . . . , s; satisfying Re(sj) + ■ • • + Re(s;) > I — j + 1 for all j = 1, . . . , I. We 
call |s| = si + ■ ■ ■ + si the weight and denote the length by l{s). The special values of multiple 
zeta functions at positive integers have significant arithmetic and algebraic meanings, whose defining 
series (QJ will be called MZV series, including the divergent ones like C(- ■ • > !)• These obviously 
generalize the notion of harmonic series whose weight is equal to 1. 

MZV series are related to many aspects of number theory. One of the most beautiful computations 
carried out by Euler is the following evaluation of zeta values at even positive integers: 

('271 "l^™ 

where are Bernoulli numbers defined by the Maclaurin series t/(e* — 1) = X^fc^o B^t^ /k\. In this 
paper, we will study partial sums of MZV series which turn out to be closely related to Bernoulli 
numbers too. These sums have been called (non-alternating) multiple harmonic sums and studied by 
theoretical physicists (see and their references). Their main interest is fast computation of their 
exact values and the algebraic relations among them. Our focus, however, is on p-divisibility of these 
sums for various primes p and special attention will be paid to the cases where the sums are divisible 
by higher powers of primes than ordinarily expected which is often related to the irregular primes, 
i.e., primes p which divide some Bernoulli numbers Bt for some positive even integers t < p — 2. We 
say in this case (p, t) is an irregular pair. 

1.1 Wolstenholme's Theorem 

Put s := (si, . . . , s/) G N' and denote the nth partial sum of MZV series by 

i/(s;n) = i/(si,...,s/;n) := ^ V'-'V. " e Z>o. (2) 

l<fei<-<fe;<n 

*2000 Mathematics Subject Classification: Primary: 11A07, 11Y40; Secondary: 11M41. 
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By convention we set H{s; r) = Q for r = 0, — 1, and i?(0; 0) = 1. To facilitate our study we 
also define (cf. ^51 OS) ) 

5(51, . . . , s;; n) := ^ k^''^ ■ ■ ■ k^"' , n e Z>o. (3) 

l<ki<-<ki<n 

To save space, for an ordered set (ei, . . . , e*) we denote by {ei, . . . , et}'^ the set formed by repeating 
(ei, . . . , et) d times and = {1}'^. For example iJ({s}'; n) is just a partial sum of the nested zeta 
value series C(s) of length I which we refer to as a homogeneous partial sum. The partial sums of 
nested harmonic series are related to Stirling numbers St{n,j) of the first kind which are defined by 
the expansion 

n 

x{x + l){x + 2) ■ ■ ■ {x + n - 1) = ^ St{n,j)x^. 

i=i 

We have St{n,n) = 1, St{n,n — 1) = n{n — l)/2, and St{n, 1) = (n — 1)!. It's also easy to see that 
St{n,j) = {n-iy.- H{l^-^;n-l), forj = I,-- - ,n. (4) 
Now we restrict n to prime numbers. 

Theorem 1.1. (Wolstcnholmc 1862 24, p. 89]) For any prime number p > 5, St{p, 2) = (mod p^). 

We find on the Internet the following generalization of the above theorem by Bruck 8 although 
no proof is given there. Denote by p{m) the parity of m which is 1 if to is odd and 2 if to is even. 

Theorem 1.2. For any prime number p > 5 and positive integer I — 1, . . . ,p ~ 3, we have 
St{p,l + l)=0, iJ(l';p-l)=0 (mod 

This of course implies that p\St{p, Z) for 1 < / < p which was known to Lagrange p. 87]. 

It is noticed that not only p"^ but also p^ possibly divides St{p, 2), though rarely, and therefore p^ 
possibly divides the numerator of iJ(l;p— 1) written in the reduced form. So far we know this happens 
only for p = 16843 and p = 2124679 among all the primes up to 12 million (see OSl [T31 1^ [TU ] ). 
The reason, Gardiner told us in |18| . is that these two primes are the only primes p in this range 
such that p divides the numerator of Bp-3. Bruck (H) further gave a heuristic argument to show that 
there should be infinitely many primes p such that p'^ divides S{p, 2), which is equivalent to say that 
there are infinitely many irregular pairs {p,p — 3). 

There are other generalizations of Wolstenholme 's Theorem to multiple harmonic sums. Bayat 
proved (corrected version, see Remark Oil 

Theorem 1.3. (4, Thm. 3]) For any positive integer s and prime number p > fc + 3 we have 

H{s;p-l)=0 (mod/'''+^^). 

In 1321 Slavutskii showed 

Theorem 1.4. f |32l Thm. 2]) Let s and m be two positive integers such that (to, 6) — 1. Let 
t — {<p{m) — l)s and v1„(to) = rip|m(l ~ P"^^)Bn be the Agoh's function, where 4>(rn) is the Euler's 
Phi function. Then 

mAt{m) (mod to^) if s is even, 

(i/2)m^j4(_i (to) (mod to^) if s is odd. 



l<k<m 
(k.p) = l 
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1.2 Homogeneous multiple harmonic sums 

The classical result of Wolstenholme is the original motivation of our study. We will prove the 
following generalization of Thm. 11.11 and Thm. 11.21 to homogeneous multiple harmonic sums (see 
Thm. imi) 

Theorem 1.5. Let s and I be two positive integers. Let p be an odd prime such that p > I + 2 and 
p — 1 divides none of Is and ks + 1 for k = 1, . . . ,1. Then 

H{{sy;p-1) =0 (mod pP^''*-!^). 

In particular, the above is always true if p > Is + 3. 

We also look at some cases when the congruences hold modulo higher powers of p. Recently, 
Zhou and Cai 39 prove that 

Theorem 1.6. Let s and I be two positive integers. Let p be a prime such that p > Is + 3 

\ ("l)' t,r4.y ^P-'^-2 (^^od/) ^f2\ls, 
H{{sy;p-l)^S{{sy;p-l)^\ 2(Zs + 2) 

[ (-l)'"'|^Sp-is-i (mod/) if2\ls. 

We will prove an analog of this in the non-homogeneous even weight length two case (see 
Thm. 1221). 

1.3 Non-homogeneous multiple harmonic sums 

The third section of this paper deals with non-homogeneous multiple harmonic sums. We consider 
the length 2 case in Thm. 13.11 whose proof relies heavily on generating functions of the Bernoulli 
polynomials and properties of Bernoulli numbers such as Claussen-von Staudt Theorem. 

Theorem 1.7. Let si,S2 be two positive integers and p be an odd prime. Let si = m,S2 = n 
(mod p — 1) where < m, n < p — 2. If m,n > 1 then 



(—1)" f rn + n 



H{si, S2;p — 1) = < TO + n \ m 



Bp-m-n (mod p) ifp>m + n, 



(mod p) ifp<m + n. 

The same idea but more complicated computation enables us to deal with the length 3 odd weight 
case completely (see Thm. I3.5|l . 

Theorem 1.8. Let p be an odd prime. Let (si, 32,53) G N'^ and < l,m,n < p ~ 2 such that 
Si = /, S2 = TO, S3 = n (mod p — 1). If l,m,n > I and w = l + m + nisan odd number then 

11(1, m,n;p ~ 1) = 1(1, to, n) — I{n, ni, I) (mod p) 

where I is defined as follows. Let w' ^ w — {p — 1) if p < w < 2p and w' — w otherwise. Then 

'0 if w > 2p, or if I + m < p and p < w < 2p ~ 1, 

/(Z,TO,n)= Jl/2n tfw^p,2p^l. 



{-ly^+^f'^')^^^ otherwise. 
. ^ ' \n) 2w' 



However, it seems to be extremely difficult to adopt the same machinery for general larger length 
cases. For the even weight cases in length 3, we are only able to determine the p-divisibility for 
iJ(4, 3, 5;p — 1), iJ(5, 3, 4;p — 1) and the three multiple harmonic sums of weight 4: if(l,l,2;p— 1), 
iJ(l, 2, l;p — 1), and H{2, 1, l;p — 1), which are distinctly different from the behavior of others. 

Recently, M. Hoffman studies the same kind of questions independently in '53', "5^^ from a different 
viewpoint. We strongly encourage the interested reader to compare his results to ours. For example, 
Hoffman defines the convolution operation on composite of indices (see §6]). We can apply this 
to a few of the above results to find more Wolstenhomles' type congruence in section 1X71 



3 



Theorem 1.9. Let p be a prime and s G N'. Assume p > |s| + 2. Then 



H{s;p~l) = S{s;p-l) = (mod p) 
provided s has one of the following forms: 

1. s ^ (l™7 2, 1") for m,n > and m + n is even. 



2. (l", 2, V'-\ 2, w/iere n>2 is even. 

3. s= 2, 2, 1") where n > 2 is even. 
^. s = (l", 2, 1", 2, 1") w/iere ti > 0. 

In j38| we will look at the mod p structure of the multiple harmonic sums for lower weights. 

One can also investigate the multiple harmonic sum H{s; n) with fixed s but varying n. We will 
carry this out in the second part of this series 36^ . Such a study for harmonic series was initiated 
systematically by Eswarathasan and Levine ^7] and Boyd [7], independently. 

The theory of Bernoulli numbers and irregular primes has a long history, and results in this di- 
rection are scattered throughout the mathematical literature for almost three hundred years starting 
with the posthumous work "Ars Conjectandi" (1713) by Jakob Bernoulli (1654-1705), see ^1]. With- 
out attempting to be complete, we only list some of the modern references at the end. In particular, I 
learned a lot from the work by Buhler, Crandall, Ernvall, Johnson, Metsankyla, Sompolski, Shokrol- 
lahi, [51 IIUI im I13L 1151 1161 128| . and Wagstaff on finding irregular primes. For earlier history, one 
can consult '35' and its references. Often in my computation I use the table for irregular primes less 
than 12 million available online at ftp://ftp.reed.edu/users/jpb maintained by Buhler. My interest 
on multiple harmonic sums was aroused by the work of Boyd 0, Bruck 0, Eswarathasan and E. 
Levine (17. , and Gardiner '18' on the nice and surprising relations between partial sums of harmonic 
series and irregular primes. I'm indebted to all of them for their efforts on improving our knowledge 
of this beautiful part of number theory. 

2 Generalizations of Wolstenholme's Theorem 

It's known to every number theorist that for every odd prime p the sum of reciprocals of 1 to p — 1 
is congruent to modulo p. However, it's a little surprising (at least for me) to know that the sum 
actually is congruent to modulo if p > 5. This remarkable theorem was proved by Wolstenholme 
in 1862. 

2.1 Generalization to zeta- value series 

To generalize Wolstenholme's Theorem we need the classical Claussen-von Staudt Theorem on 
Bernoulli numbers (see, for example, ^3 p. 233, Thm. 3]): 

Lemma 2.1. For m G N, -62™ + X]p-i|2m ^/P ^""^ integer. 

We begin with a special case of our generalization which only deals with zeta- value series, i.e., 
MZV series of length 1. The general case will be built upon this. 

Lemma 2.2. Let p he an odd prime and s he a positive integer. Then 



Remarks 2.3. (1) The conditions in Bayat's generalization of Wolstenholme's Theorem ^ Thm. 3] 
should be corrected. For example, taking k = 2 and p = 5 in 4, Thm. 3(i)] we only get iJ(3;4) = 
2035/1728 ^ (mod 5^). In general, if 2fc =p-l in 4, Thm. 3(i)] we find that H{2k-l]p-l) = -p 
(mod p^) by taking s — 2k — 1 and n = 1 in our lemma. 

(2) Most of the lemma follows from Slavutskii's result Thm. ^31 However, we need a more direct 
proof which we will reference later. 




(modpP("+i)) ifp-l\s,s + l, 

(mod p ) if p — l|s, 

(mod p^) if s + 1 ~ n{p — 1). 
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Proof, li p — l\s then H{s;p — 1) = p — 1 = — 1 (mod p) by Fermat's Little Theorem. So we assume 
p — 1 \ s. This implies that the map a — > a^'' is nontrivial on (Z/pZ)^, say b" ^ 1 (mod p)for 
some 1 < < p . Then it's not hard to see that (1 — b^'')H{s;p — 1) = (mod p) and therefore 
H{s;p — 1) = (mod p). This holds for any s such that p — 1 \ s, whether it is even or odd. 

The last case of modulus p^ for odd s can be handled by the same argument as in the proof of 
Wolstenholme's Theorem. We produce two proofs below for both completeness and later reference. 

Let s be an odd positive integer. Choose n large enough so that t := np{p — 1) — s > 3 is odd. 
Then by the general form of Fermat's Little Theorem 

k=l k=l k=l 

By a classical result of sums of powers (see [23 P- 229]) we know that 

P-i ^ 

Y^k* ^-—{Bt+i{p)-Bt+i) fort>l, (5) 
fc=i 

where Bm{x) are the Bernoulli polynomials. Further, 

b,^^[p)JY^('^Ab,p^+'-^. 

j=o \ / 

Observing that pBj is always p-integral by Lemma |2. II we have 

Y,k' =pBt + -p^Bt^i (modp2). (6) 

k=l 

When p — 1 \ s + 1 the lemma follows from the facts that Bj = if j > 2 is odd and that 
Bt-i = Bnp(p-i)-s-i is p-integral by Lemma I^TTl If p — l|s + 1 and s + 1 = m{p — 1) then we choose 
n = m in the above argument. Then t = mp{p — 1) — m(p — 1) + 1 is odd. So Bt = 0, p — l\t — 1 and 
pBt-i = — 1 (mod p) by Lemma l2 . II again. From congruence 10) we get 

H{s]p - 1) = -pt/2 = -p{n + l)/2 (mod p^). 

In fact, there's a shorter proof for the odd case which is not as transparent as the above proof. 
By binomial expansion we see that 



A, — X f\, — X l\, — X 



Noticing that l/{p — ky = — l/fc" (mod p) we have 

fc=i k=i 

whenever p — 1 \ s + 1 because p divides X]fe=i l/fc"^^ by the even case. □ 

When s = p*^ we can work more carefully with binomial expansion in the shorter proof and see 
that p2+e cjivides H{s;p— 1). When e = 1 this explains the fact that 125 divides H{5; 4). We record 
the phenomenon in the following proposition. 

Proposition 2.4. Let s be a positive integer and Vp{s) = v and Vp{s + 1) = u (so that either u or v 
is 0). If p > 5 is a regular prime then 



Vp{H{s;p^l)) 




ifp- l\s, 

if s is even and p — 1 \ s, 
if s is odd and p — 1 |.s + 1, 
if s is odd and p — 1 \ s + 1. 
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Suppose p is irregular and let 1 < m < p{p — 1) such that 

{— s (mod p{p — 1)) if s is even, 

—s — 1 (mod p{p — 1)) if s is odd. 

If p — l|s + 2, s + 3 and p^ \ Bm/ni then the nonzero valuations can increase by at most 1. 

Proof. We only consider the case when p — 1 \ s. Suppose s — p'"a and p \ a. Let e > u + 2 be any 
positive integer such that t := p'^{p — 1) — s > 1. Then p ~ 1 \ t. It foUows from Fermat's Little 
Theorem that 

p-i 

iJ(s;p- 1) EE Vfc* (modp''+3). 



fe=i 



If s is even then t is even and 



H{s- p-l) =pBt + ^-^^-^p^Bt^2 (mod 
6 

_ v+l^t' (^-l)(i-2) „+3 , „+3, 

= -ap ^ — H — p ^ Bt" (mod p ^ ) 

by Kummer congruences, where t = t' ,t — 2 = t" (mod p{p — 1)) and 1 < t', t" < p{p — 1). Note that 
Wp-valuation of the first term is w + 1 or higher depending on whether {p,t') is regular pair or not. 
The smallest Up-valuation of the second term is w + 2 which happens if and only if p — If p is 
irregular then p — I \ t" hy our assumption p — I \ s + 2 and hence the second term is always divisible 
by p""*"^ . It follows that the tip-valuation is w + 1 if p is regular and is at most w + 2 if p is irregular 
since we assumed p^ | Bfi /t'. This proves the proposition when s is even. 

When s is odd, then we need to consider two cases: u = or v = 0. Both proofs in these two 
cases are similar to the even case and hence we leave the details to the interested readers. □ 

Remark 2.5. One can improve the above by a case by case analysis modulo higher p-powers. For 
example, one should be able to prove that the nonzero valuations can increase by at most 6 if p is 
irregular less than 12 million. 

Remark 2.6. Numerical evidence shows that if p is irregular then the nonzero valuations seems to 
increase by at most 1 in most cases. The first counterexample appears with p — 37, and s — 1048. 
Note that s is even and p— 1 f s so that if one leaves out the conditions in the lemma then the prediction 
would say Vp{H{s;p — 1)) is at most 2 because 37 is irregular. But we have U37(_ff (1048; 36)) = 3 
because 37^|i?p(p_i)_s = -B284- 

Corollary 2.7. Let s > A be a positive integer. Let p > 3 be a prime. If p is irregular then we 
assume it satisfies the conditions in the proceeding proposition. Then H(s;p — 1) ^ (mod p"). 

Proof. The case p = 3 can be proved directly because 1 + 2* < 3^* if s > 2 and 

1 + — = ^— ^ (mod 3^). 

Suppose p > 5. Let s — p" a where p \ a. If p is regular then by Prop. 12.41 the largest value of 
Vp{II(s;p — 1) is w + 2 which is less than s because s > 4 and p > 5. If p is irregular satisfying the 
conditions in Prop. then the largest value of Vp{II{s;p — 1) is at most v + 3 which is still less than 
s because s > 4 and p > 37. □ 



2.2 p-divisibility, Bernoulli numbers, and irregular primes 

Numerical evidence shows that congruences in Lemma 12.21 is not always optimal. For any zeta value 
series, every once in a while, a higher than expected power of p divides its (p — l)-st partial sum. A 
closer look of this phenomenon reveals that all such primes are irregular primes. Going through the 
proof of Lemma 12.21 a bit more carefully one can obtain the following improvement. 
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Theorem 2.8. Suppose n is a positive integer and p is an odd prime such that p > 2n + 3. Then 
we have the congruences: 

—2 2n 

H{2n-l-p-l)=p-H{2n-p-l)=p'' ■——■Bp^2n-i (mod p^). (7) 



2n-l ' ' ' ' 2n + l 

Therefore the following statements are equivalent: 

(1) . Sp_2„-i = (mod p). 

p-i 

(2) . i/(2n;p- 1) =^l/fc2" = (modp2). 

fc=i 

(3) . i/(2n- l;p- 1) = = (mod p^). 

A;=l 

(4) . H{n,n;p-l)^ ^ l/(fcifc2)" = (mod p^). 

l<fel<fe2<P 

Proof. The congruence relation Q and the equivalence of (1) to (3) in the theorem follows from the 
two congruences in [301 after |E1 (16)]- See also ^01 P- 281]. The equivalence of (2) and (4) follows 
immediately from the shuffle relation*: 

H{n;p- 1)2 = 2H{n, n;p-l) + H{2n; p - 1). (8) 

□ 

Remark 2.9. This theorem also follows directly from a result of See Thm. 

Gardiner ^B] proves the special case of the equivalence when n — 1. He has one more equivalence 
condition which involves the combinatorial number (^) . Other variations of the classical Wolsten- 
holme's Theorem can be found in 1 . It would be an interesting problem to find the analog for the 
zeta value series. It is also worth mentioning that it is not known whether (^^^ = 2 (mod n'^) would 
imply n is a prime, which is called the converse of Wolstenholm's Theorem |31| . 



2.3 Hoffman's convolution 

In this section we recall Hoffman's convolution operation of the composites and provide new con- 
gruence modulo a prime square. These results will be extremely useful when we deal with multiple 
harmonic sums of arbitrary lengths. 

Let's first recall the definitions. Let A: be a positive integer and s = (ii, . . . , i^) of weight n = |s|. 
We define the power set to be the partial sum sequence of s: P{s) = {ii,ii + i2, ■ ■ ■ ,ii + ■ ■ ■ + ik-i) 
as a subset of (1, 2, . . . , n). Clearly P provides a one-to-one correspondence between the composites 
of weight n and the subsets of (1, 2, . . . , n — 1). Then s* is the composite of weight n corresponding 
to the complimentary subset of P{s) in (1, 2, . . . , n — 1). Namely, 

s* =P-I((l,2,...,n-1)-P(s)). 

It's easy to see that s** = s so * is indeed a convolution. For example, if ii^ik > 1, «2i ■ ■ • j *fc-i > 2 
then we have 

(zi, . . . , ik)* = (i'^"\ 2, 2, . . . , 2, V"-'-^, 2, V"-'). 

Further we set the reversal s = (s;, . . . , si). The following important result is due to Hoffman: 
Theorem 2.10. f |25[ Thm. 6.8]) For s ~ (si, . . . , s;) and any positive integer n define 

S{s-n)^ J2 (9) 

l<fel<-<fci<Tl 

*Somc authors call this a stuffle relation, i.e., "shuffling" plus "stuffing". 
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Then for all prime p 

Sis*;p-1)= -S{s;p-1) (modp), (10) 

5(s;p- 1) =(-l)l"IS'(s;p- 1) (modp). (11) 
We also have the following equalities: 

5(s;n)=^iJ(i;n), (12) 

i^s 

where i ^ s means i can be obtained from s by combining some of its parts, and 

I 

(-l)'(^)5(s;p-l)= i-^Yli^^^J'P-^) (13) 

where LJj"=i catenation of Si to S;. 

Proof. Note that in 25^ the right hand side of © is denoted by S{s; n) instead. But it's not difhcult 
to verify that for any s 

s* = s^. (14) 

So (|10|l is equivalent to 25, Thm. 6.8]. Congruence pi|) foUows readily from the substitution of 
indices: k p — k. 

Equations and p2|) relate S'-version multiple harmonic series and our TJ-version (denoted 
by A by Hoffman). Equation (|12|) is Eq. (7)]) and is equivalent to the second unlabelled 
formula in the proof of 25, Thm. 6.8]. □ 

Next we want to provide a more precise version of congruence (jlUII . 
Theorem 2.11. Let s be any composite of weight w. Let p be an arbitrary odd prime. Then 

-S{s*;p-l)=S{s;p-l)+p^H{tU{l};p-l)^ (mod p^). (15) 
Proof. For any sequence {f{n)}n>-i, f{~^) = know there are two operators S and V: 

71 

S/(n) = 5]/(^), Vfin) = f{n)-fin-l), Vn > 0. 

i=0 

Recall that 'SVS{s; n) = —S{s*; n) for any s and positive number n by [261 Thm. 4.2]. Then for any 
odd prime p we have 

P 

U: 

=0 
/p-2 



-S{s*:p- 1) =SV5(s;p - l){n) - ^ ( J i j (-^TSi^' 



e5(s; p-l){n)+p z)^ (mod 

.S{s-p - l)(n) + p I g E I (-od P^) 

E5(s;p-l)(n)+p K]i/(tU{l};p-l) (modp^). 



, t-<s 



□ 
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2.4 Wolstenholme type theorem for homogeneous sums 

In the above we have studied the p-divisibility of H{s;p — 1) for positive integers s. It is difficult to 
start with arbitrary s so we only consider the homogeneous multiple harmonic sums at the moment. 
First we can easily verify the following shuffle relations: for any positive integers n, m, si, . . . , s;, 

I 

H{m;n) ■ H{si, . . . , s;;n) = ^ H{s;n) + ^iJ(si, . . . , Sj_i, Sj + m,Sj+i, ...,si;n) 

seShfl(m,(si, •••,«;)) J = l 

where for any two ordered sets (ri, . . . , rt) and {rt+i, . . . ,rn) the shuffle operation is defined by 



Shfl((ri,...,rf),(rt+i,...,r„)) := |J ('■<t(i), 



tr permutes {l,...,n}, 
a-\l)<---<<7-\t), 
o--^(t+l)<---<<T~i(n) 



Hence for any fc = 1, . . . , Z — 1, we have 



H{ks-n) ■ H{{sy-'';n) ^ ^ H{s;n) + ^ H{s;n). 

seShfi ({fcs},{s}'-'=) seShfl ({(fe+i)s},{s}'-'«-i) 



Applying X]a;=\(^1)'^ ^ on both sides we get 
Lemma 2.12. Let s,l and n he three positive integers. Then 

I 



lH{{sy;n) ^Y.^-lf'^H(ks-n) ■ H{{sy-^-n). (16) 
fc=i 

We now need a formula of the homogeneous multiple harmonic sum iJ ({s}'; n) in terms of partial 
sums of ordinary zeta value series. Let P{1) be the set of unordered partitions of /. For example, 
P(2) = {(1, 1), (2)}, P(3) = {(1, 1, 1), (1, 2), (3)} and so on. 

Lemma 2.13. Let s,l and n be positive integers. For X — (Ai, . . . , A,.) G P{1) we put H\{s;n) = 
Yii^i H{XiS;n) and define c\ = 1)*^ ' n[=i(-^« ~ Then 

nH{{sy;n) = J2 CAi?A(s;n). (17) 
AeP(0 

Proof. It follows easily from ^26, Thm. 2.3]. □ 

It is obvious that C(i ... i) = 1. When I = 3 we have C(i^2) — — 2ci + C2 — —3, which implies 

6iJ ({s}3; n) ^ H{s- nf ~ 3H{s; n)H{2s; n) + 2H{3s; n). (18) 

When I = 4: we have C(i_i^2) = C(i_2) - 3c(i_i) = -6, C(i 3) = C3 + 6ci = 8, and C(2,2) = 3c2 = -3, 
which implies 

24:H{{s}'^; n) ^ H{s; n)^ - 6H{s; nfH{2s\ n) + 8i7(s; n)H{3s; n) - 3H{2s; nf - QH{As\ n). 

Theorem 2.14. Let s and I be two positive integers. Let p be an odd prime such that p > I + 2 and 
p — 1 divides none of si and ks + 1 for k ~ 1, . . . ,1. Then the homogeneous multiple harmonic sum 

H{{sy;p- 1) = 5({s}';p - 1) = (mod pP^'^-^^). 

In particular, if p > Is + 3 then the above is always true and so p|iJ({s}';p — 1). 

Proof. Congruence for H follows from Eq. H17|) and Lemma Congruence for S then follows from 
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Recently, Zhou and Cai obtain an improved version of the above result, see |39) 
Theorem 2.15. Let s and I be two positive integers. Let p be a prime such that p > Is + 3 

H{{sy;p~ 1) ^ S{{sy;p - 1) ^ \ 2(^5 + 2) 

[ (-l)'-'lj^Sp-is-i (modp2) ,f2\ls. 

Proof. Congruence for H follows from |39| . Congruence for S then follows from and an induction 
on I. □ 

2.5 Higher divisibility and distribution of irregular pairs 
From Thm. ITT^ we immediately have 

Proposition 2.16. Let p be an odd prime. Suppose {p,p — Is — 2) is an irregular pair (so both I and 
s must be odd numbers). Then 

H{{sy;p- 1) = S{{sy;p- 1) = (mod p'). 

Suppose {p,p — Is — 1) is an irregular pair. Then 

H{{sy;p~ 1) = S{{sy;p- 1) = (mod p'). 

One of the first instances of this proposition is when s ~ 1 and I ~ 3 given by the first irregular 
pair (37,32), which, by our proposition, implies that _ff(l,l,l;36) = (mod 37'^). The case when 
s = 3 and / = 3 appears first with the irregular pair (9311,9300). So we know iJ(3, 3, 3; 9310) = 
(mod 9311^). 

We believe these are just the first two of infinitely many such pairs because evidently not only 
Bernoulli numbers but also the difference p — t for irregular pairs {p, t) are evenly distributed modulo 
any prime. Precisely, we have the following 

Conjecture 2.17. For any fixed positive integer M and integer c such that Q < c < M we have 

iUp.t) : prime plBfA even,p < X,n ^ t = c (mod A/)| | , / I ^ I ' 

lim ' ^ ^[ '-f- — ^ = I I M if2\M, (19) 

x~,oo i{{p,t) : prime p Bt,t even,p < X} 

[2/A/ */2f c, 2\M. 

Further we can replace the sets by restricting p to all irregular primes with a fixed irregular index 
which is defined as the number of such pairs for a fixed p. 

As a result we expect there are about one-third irregular pairs satisfying the conditions in 
Prop. [TTEl In Table^we count the first 11,000 irregular pairs. We denote by N{k,m) the number 
of irregular pairs {p, t) satisfying p — t = k (mod 3) in the top m irregular pairs, < fc < 2, and 
by P{k,m)% the percentage of such pairs. For irregular primes with fixed index we compiled some 
more tables in the Appendix at the end of this paper. 

3 Non-homogeneous sums 

Having dealt with homogeneous multiple harmonic sums we would like to do some initial experiments 
on the non-homogeneous ones. 
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m 


1,000 


2,000 


3,000 


4,000 


5,000 


6,000 


7,000 


8,000 


9,000 


10,000 


11,000 


N{0,m) 


322 


664 


996 


1318 


1637 


1978 


2310 


2628 


2968 


3310 


3672 


P(0,m) 


32.20 


33.20 


33.20 


32.95 


32.74 


32.97 


33.00 


32.85 


32.98 


33.10 


33.38 


N{l,m) 


343 


666 


1005 


1352 


1676 


1999 


2340 


2683 


3017 


3354 


3673 


P{l,m) 


34.30 


33.30 


33.50 


33.80 


33.52 


33.32 


33.43 


33.54 


33.52 


33.54 


33.39 


7V(2,m) 


335 


670 


999 


1330 


1687 


2023 


2350 


2689 


3015 


3336 


3655 


P(2,m) 


33.50 


33.50 


33.30 


33.25 


33.74 


33.72 


33.57 


33.61 


33.50 


33.36 


33.22 



Table 1: Distribution of p — t (mod 3) for irregular pairs {p, t). 



3.1 Non-homogeneous sums of length 2 

We begin with C(l, 2) and <^{2, 1) scries. For any positive integer n from the shufHe relation we have 

H{1, 2; n) + H{2, 1; n) = H{1; n) ■ H{2; n) - H{3; n). (20) 

However, it seems to be quite difHcnlt to disentangle H{1, 2; n) from H{2, 1; n). Maple Computation 
for primes up to 20000 confirms the following 

Theorem 3.1. Let si,S2 be two positive integers and p be an odd prime. Let si = m,S2 = n 
(mod p—l) where < m,n < p — 2. Then 



H{si,S2;p- 1) = < 



1 



(-1)" fm + ri 



(modp) z/(m,n) = (0,0), (1,0), 
(mod p) if (to, n) = (0, 1), 



TO + n V TO 




B. 



p—m—n 



(mod p) if p > m + n and m,n> 1, 
(mod p) otherwise. 
In particular, ifp>s + t for two positive integers s and t then we always have 

i-iy fs + f^ 



S{s,t;p-l) = H{s,t;p-l) 



s + t 



B„ 



-t (modp). 



Proof. We leave the trivial cases to the interested readers and assume in the rest of the proof that 
TO, n > 1. 

Let M = p - 1 - m and N = p-l-n. Then 1 < M,N <p - 2. By Fermat's Little Theorem 

H{m,n;p-1)= ^ /cffc^ (mod p). 

l<fel</S2<p-l 

Define the formal power series in two variables 



r,s=0 \ l</ci<fe2<f)-l 



r\s\ 



Exchanging summation we get 



E 



f{x,y)^ > e 

1</Sl<fc2<p-1 

oo oo 



kix+k2y 



pp{x+y) _ pX+y 



{ePy - e2')e^ 



(e^+2' - l)(e^ - 1) {ev - l){e'= - 1) 



B,+,{p) - B,+,{1) 



ip(B,+ix^(x + y)^-B,+i(l)xV 



where B„(l) = if to 7^ 1 and Bi{l) = -Bi = i. We only care about the above sum when 
1 < p — 2 since 1 < M, N < p — 2. So we may as well replace 5^(1) by B^ everywhere in the 
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last display and throw away the terms with j = 0, ±1. The resulting power series is 

j=2 1=1 ^ \ / 

Let I + i — M and j — I = N. Then we see that the coefficient of x'^'^y^ is 



M 



E 



-^{M -l + l)\{N + I 



N + l 



E 



pBm-i+iBn+i 



^ (M - z + ly.miv 



(21) 



Note that 0<M-; + l<M<p-2and2<iV+l<Af + ;<M + iV<2p-4. Ifi\f + Ar<p-2 
(i.e. m + n > p) then H21|l is always congruent to mod p. Otherwise, if M + > P — 2 then all the 
terms are in pZp except when N + I = p ^ 1. Hence 



MINI 



pB 



M+N+2-pBp-i 



l<ki<k2<p-l 

So finally we arrive at 



{M + N + 2 - py.N\{p - 1 - TV)! 



H{'m, n;p — 1) 



E 



l<fcl</C2<p-l 



k™k~ 



1 '^2 



-{p-m- ly.Bp^yn^r 

n\(j) — m — ny. 



(mod p). 



(mod p). 



One can now use Wilson's Theorem to get the final congruence for H in our theorem without too 
much difficulty. For 5 we now use S{s, t) = H{s, t) + H{s + t) and Wolstenholme's theorem. □ 

Taking m — l,?i = 2 in the theorem we obtain H{l,2;p — 1) = i?p_3 (mod p) for p > 3. We 
verified this on Maple for the only two known irregular pairs of the form (p,p — 3), namely, p ~ 16843 
and p = 2124679. If we take m = 2, n = 3 we find that H{2, 3;p — 1) = i?p-5 (mod p) ii p > 5. 
There is only one irregular pair of the form {p,p ~ 5) among all primes less than 12 million, namely, 
(37, 32). Indeed, for p = 37 computation shows that H{2, 3; 36) = (mod 37). As a matter of fact, 
we formulated Thm. mi onlv after we had found these intriguing examples. 

We now can provide an analog of Thm. lO^ in the non- homogeneous case of length 2. 

Theorem 3.2. Let p be an odd prime. Suppose s and t are two positive integers of same parity such 
that p > s + t + 1. Then 



H{s,t;p- 1) ^p 
S{s,t;p- 1) =p 



{-lyt 

i-iyt 



s + t + l 
s 

s + t + l 
s 



-i-iys 
-{-lys 



s + t + l 

t 

s + t+l 

t 



- s-t 
+ s + t 



Bp-s-t-1 
2{s + t+l) 



Bp-s-t-i 
2{s + t+l) 



(mod p^) 
(mod p^). 



Proof. By the shuffie relation (dropping p ~ 1 again) we see that 

H{s) ■ H{t) ^ H{s, t) + H{t, s) + H{s + t). 
By the conditions on Si and S2 we know from jT)) 

H{s) ■ H{t) EE 0, H{s + t)= ^^'^ ^ ^\ B 



Therefore 



H{s,t)+H{t,s) 



s + t + l 
-p{s + 1) 



s + t + l 



p-s-t-i (mod p ). 
Bp-s-t-i (mod p^). 



(22) 

(23) 
(24) 
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Moreover, by the old substitution trick i,j p — i,p — j 

1 



= E + + -)* (rnodp^) (s + tiseven) 

l<j<i<p * * 



i<j<j<p 



Ei? (t, s) + psH{t, s + 1) + + 1, s) (mod p^ 



=H{t,s)+p 



s + t + 1 



s + t+l 

t + 1 



Bp-s-t-1 
s + t + 1 



(mod p^). 



Combined with (|24|l this completes the proof of the congruence for H. Then the 5* part follows from 
the identity S{s, t) = H{s, t) + H{s + t). □ 



3.2 The case of palindrome s 

Recall that for s = (si, . . . , s;) we have set its reversal s — (s/, . . . , si). 

Lemma 3.3. Let p be an odd prime. Let I be a positive integer and s G N'. Let |s| = X]i=i 
weight of s . Then 

H{s;p- 1) =(-l)l^li/(s;p- 1) (mod p), 
S'(s;p- 1) =(-l)l"l5(s;p- 1) (mod p). 

Proof. Use the old substitution trick ki ^ p — ki for all i in the definitions Q and ||2J). □ 

An immediate consequence of this lemma is 

Corollary 3.4. Let p be an odd prime. Lf s = s and |s| is odd then 

H{s;p-l) = S{s:p-l) = (mod p). 

On the contrary, a lot of examples show that if the weight |s| > 6 is even and if the length is 
bigger than 2, then we often have H{s;p — 1) ^ (mod p) when p is large (say, p > 2\s\), even in 
the case that s is a palindrome. For example, in length 3 if s 7^ (4,3,5), (5,3,4) then this seems to 
be always the case (see Problem 13. in|l . A remarkable different pattern occurs for length 3 weight 4 
case which we will consider in subsection 13.41 Clarifying this completely might be a crucial step to 
understand the structure of i7(s;p — 1) in general. 



3.3 Multiple harmonic sums of length 3 with odd weight 

One may want to generalize Thm. IXTl to multiple harmonic sums of longer lengths. However, the 
proofs become much more involved. Extensive computation confirms the following 

Theorem 3.5. Let p be an odd prime. Let (si, 32,53) G N'^ and < l,m,n < p — 2 such that 
Si = I, S2 = rn, S3 = n (mod p — 1). Then 

H{0, m,n\p — \) = H{in — 1, rt;p — 1) — H{m, n;p — 1), (mod p), 

H{l,0,n;p-1) = H{l,n-l;p-l) -H{l,n;p'- 1) H{1 - l,?i;p- 1) (modp), 
H{1, m, 0;p — 1) = ~H{1, m — l;p — 1) — H{1, m;p — 1) (mod p). 

Suppose l,m,n> 1. Suppose further that w = l + m + nisan odd number. Then 

H(l, m, rt;p — 1) = 1(1, m, n) — I{n, m, I) (mod p) 



13 



where I is defined as follows. Let w' ^ w — {p — 1) if p < w < 2p and w' = w otherwise. Then 
I{l,m,n) = / ,\ rj 




if w > 2p, or if I + m < p and p < w < 2p — I, 
ifw=p,2p- 1, 



\n J 2w' 

In particular, if a prime p > I + m + n for positive integers I, m, n such that w = I + m + n is odd 
then 

^ / w \ I w ^ 
H{l,m,n;p ~ 1) = —S{l,m,n]p— 1) = (—1)" 



I J \n 

Remark 3.6. In j26[ Thm. 6.2] Hoffman gives a easier proof tlian ours 



™ (mod p) . 



2w 



Proof. Tlie last congruence follows easily from relations Hll|l and H13|l so we only consider the first 
part. Throughout the proof all congruences are modulo p. We leave the trivial cases to the interested 
readers and assume in the rest of the proof that l,m,n > 1 and w — l + m-\-n is odd. Similar to the 
proof of Thm. Owe put L = p-l-l, M = p-l-m, and N = p-l-n. Then 1 < L,M,N <p-2. 
By Fermat's Little Theorem 

Hisi,S2,S3;p-l) = Hil,m,n;p^l)= ^ k^k^'k^. 

I<fel<fe2<fc3<p-1 

Define the formal power series in three variables 

oo I \ r s t 

r,s,t=0 yi<ki<k2<k3<p-l J I<fci<fc2<fc3<p-1 

gP(x+y+z) _ ^x+y+z ^x+y^^pz _ ^z^ 



(e=^ - l){e'^+y - l)(e^+!'+^ - 1) (e^ - l)(e=^+2' - l)(e^ - 1) 

(gPto+z) _ ^y+z-^ ^x^y ^^pz _ 

~ (e^ - l){ey - l){ey+^ - 1) ^ (e^ - l){ey - l)(e^ - 1) 
Bi+iBj+i{l)iBk+iip) - Bk+ijl)) 



oo oo 



, V V B^+i{l)Bj+,{l){Bk+iip) - B,+,{1)) , ^. , 

{^ + mJ + mk + ly■ ^ ■ 

We are interested in the x^y^'^ z^-term of the above sums where L,AI,N > 1. With this in mind 
we can safely replace all the -Bm(l) by plain B„i. It follows that H{l,m,n;p — 1) (mod p) is the 
coefficient of x^y^^z^-term of the function LlAIlNlg{x,y, z) where 



oc oo 



As i + A/ + = 3(p — 1) — w is odd, we know that for a term in (|25|l to provide a nontrivial 
contribution to x^y^ it is necessary that either j = 0, or j = 0, or fc = 1 because the only nonzero 
Bernoulli number with odd index is Bi = —1/2. But fc = 1 is not an option because all terms 
corresponding to x°-y^z cancel out. So we're left with only two cases (I) z = and (J) j = 0. We 
now handle them separately. 
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(I) When i = Q only the following terms really matter: 

^ pBiBk (x + y + z)'' - ^ pBiBj+iBk , r ^ ^ .k ki 

y — r, ^ / —r- — ,M, t + yj \[x + y + z) - z \ 

fc>2 ^ j>l,fe>2 ' 

V- yB^Bk-x'^y^-^-'^z'^-' ^ ^ pB^B^+iBk ■ x'^+^yJ+b-c^k^g-b 

^ r ■ a\(r - 1 - a)\(k - r)\ 2^ 2^2^ (j + l)c\(j~c)\a\b\ik-a-by. ^ ' 

where the first sum of (|26|l comes from setting j = —1. 

(/o) In the first sum of H2t)|) putting L — a, M = r — 1 ~ a, and N = k — r. Then we obtain the 
unique term (after multiplying L\AI\N\) 

pBiBl+m+n+1 _ pBiB^p_i)_^+i 
L + M + 1 ^ 2p-l-w + n ' 

It is an easy matter to see that the denominator is a p-unit unless I + m = p — 1 and the numerator 
is in p7jp except when w = 2p — 1 ot w = p. So we find the contribution to the coefficient of 
x^y^ z^-term in L\MlN\g{x,y, z) from this case: 



pi?li33(p_l)_u,+l 

-'0 ~^ — 'i ; — ^ ^ 

Zp — 1 — w + n 



ifu;^p,2p-l, 



if w = 2p — 1 , 
if w = p, 



p — 1 + n 

BiB2p-i-n ifl + m=p-l. 

Now let's turn to the second sum in l|26(l . Setting L = a + c, M = j + b ~ c, and N = k — a 
we see that the contribution to the coefficient of a;^y*^z^-tcrm is 

J, pBiBL+M+i-rBN+rLlMl 



(L + M + 1 - r) • (L - a)!a!(Af - ' 



Simple combinatorial argument shows that 



L\Ml ^ /A fM\ fL + M 



So we get 



^ (i - a)!a!(M - ^ \a \b 

a+b=r ^ ' ^ ' a+b=r ^ ^ 



jl _ ^t^^^ pBiBi^+M+l-rBN+ r [L + M 

2-^ L + M +l-r 

r—l 



r 



li N + L + M + 1 < p — 2, i.e., w > 2p, then all terms in /' are in pZp and we have /' = 
(mod p). Suppose w < 2p. It's obvious that N + r < L + M + N + 1< 3{p ~ 1) and L + M+l-r < 
2p — 1 — m — n < 2p — 2 . Consequently all terms in /' are in pZp except when (/i) L + M+l — r = p—1, 
or (12) N + r = p — 1, or N + r — 2{p — 1). Note that Ii and I2 can occur at the same time if 
and only if w = p. 

(/i) Suppose w < 2p and L + M + 1 — r = p — 1, i.e., r=p — I — m> 1. This occurs if and only 
if the foUowings are satisfied: I + m < p — 1 and N + r — 2p — 2~w>l (so automatically w < 2p). 
Then the corresponding term in / is 

r - RR f^ + ^^\-RR f^p-^-l-m 

However, /i = unless I + m = p — I or w = p. Hence 

-BiBp-n iil + m = p— 1 and n > 2, 

J _ J —BiBp^i ~ Bi \il + m — p— \ and 71 = 1, 

■r r- if ti; = p and n>2. 

nyn - 1) 
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(I2) Suppose w < 2p and N + r = p — 1. Because L + M > r we see that w < 2p — 2. Then the 
contribution is 



-BiBl+m+n+2-p ( L + M 



L + M + N + 2-p\p~l~N 
see case (/i) 




p — n 
(-1) 



j^^j^ BjJhp-i-w fw + 1 



2p — \ — w 



-BiB2p^i^w f2p - 2 - I ~ m 
2p — 1 — w \ n 

if w = p, 

if Z + m < p — 1 . or if p < ui < 2p, 
if / + m = p — 1 and n > 2, 

if Z + m > p and p < w < 2p — 1, or if w < p. 



(/a) Suppose w < 2p and N + r = 2(p— 1). This occurs if and only if w < p. Then the the 
corresponding term in / is 



-BiB 



L+M+N+2' 



L + M + N + 2-P 



L + M 
r 



p — w \n ~ 1 



by Kummer congruence. 

Putting Jo, Ii, I2 and Iz together we get 



(-1) 

l/2n 
(-1)"+! 
2^ 



n J2{w + 1) 



B 



p—w 



if u) > 2p, 

if p < w < 2p — 1 and I + m < p, 
if Z + m > p and w < 2p — 1, 
if w = p or w = 2p ~ 1, 
if w < p. 



(27) 



We notice the miracle that when I + m = p — 1 and n > 2, all the contributions from Iq, h and I2 
cancel out because of Kummer congruences! 

(J) When j = in (|25|l only the following terms really matter: 



E 



fe>2 



pBiBk {x + y + zf - jy + z)'- 
fc! X 



EpBlBi^lBk j / , , 



i>l.fe>2 



(i + l)!fc! 



E 



1^,6 ^k — a — b 



l<b<a<k 



pBiEkX^'-'y^z 
a\b\{k-a-b)\ 



E E 

j>l,k>2 a+b=l,a>l 



pBiBi+lBk 



{i + iy.a\bl{k - a - by. 



(28) 



(Jo) In the first sum of ^ put L = a - 1, Af = 6, and fc - a - b. After multiplying L\M\N\ 
we find the coefficient of x^y^^ z^ is 



Jo:= 



pBiB; 



3(p-l)-iii + l 



p-Z 





-1 

"2r 



if w 7^ p, 2p — 1, 
if w = p, 2p — 1. 



Let's turn to second sum in (|28|l . Setting L — i + a, M — b, and N — k — a — bwe see that the 
contribution from this sum is 



a>l 



(L + l-a)!a! 



If L + Af + + 1 < p — 2, i.e., w > 2p, then all terms arc in pZp and therefore J' = (mod p). 
Suppose w < 2p. Clearly L + l- a<p-2 and M + N + a < AI + N + L + 1< 3{p - I). We know 
that all terms in J' are in pZp except when ( Ji) Af + iV + a = p — 1, or (J2) AJ + + a = 2(p — 1). 
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(Ji) Suppose w < 2p and M + N + a = p ~ 1. This can occur if and only if M + iV < p — 2, 
i.e., m + n > p, and L + 1 — a > 1 which is equivalent to w < 2p — 2. If this is the case, then the 
contribution is 



Ji 



-BiB2p-i-^,{p-l-l)\ 
{2p-\- w)\{m + n-p+l)\ 



(-1) 



iBiB2p-i-w (w + 1 



2p — 1 — w 



I 



(J2) Suppose w < 2p and M + N + a = 2{p — 1). This occurs if and only if i + 1 — a = 
L + M + N + 3 ~ 2p ^ p — w > 1, i.e., w < p. In this case a = m + n and the corresponding term in 
B is 

_ -BiBp^^{p -l-iy. _ ^ ^y BiBp- yj ( 



^ (p — w)\{m + n)\ 
Combining cases (Jo); (^1) and (J2) we have 

'0 




p — w \l 



(-1/ 
-1 

^ , 
(-1) 



w + 1\ B2P-1- 



l J2{w + l) 



Br, 



\, 2w \l 

The theorem now follows from an easy simplification process 



if w > 2p, 

ii p < w <2p and m + n < p, 
if w < 2p — 1 and m + n > p, 

if w = p,2p — 1, 

if w < p. 



□ 



In the case when w = r + s + t is even the expression is much more complicated. For future 
reference we provide the following computation when the prime p > w + 1. The method is similar 
to that used by Hoffman to compute the length 2 case in |^. Recall that by Bernoulli polynomials 
we have 



1=1 a=0 

So modulo p we have by Fermat's Little Theorem 



d+l\ B, 



(29) 



- 1) ^ ^ 1 ^/-i- 5: p-i- 
i=i j=i 

p—l ^ i—l p—1 — 

~ E 'jt E E 

1=1 j=l a=a 
p-1 p-l-r , \ D 



P Ba jK,(a)+p-r-s~a 



a J p — r 



1=1 a=0 
p—l p—l—r 

E E 

K{a)-\-p-\-l — r — s — a 

E 

6=0 



a J p ^ r 



)-\-p—r—s — a 



a J p — r 



K{a) + P+1— r — s — a 
b 



K(a) + P+1 — r — s — a' 



where K(a) = Oifa<P~'"^s and K{a) = p — 1 if a > p — r — s. Now take the sum of powers of i 
first and observe that = (mod p) unless I = (mod p—l). Note that it is impossible to 
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have K{a) + p+ l — w — a — b = (mod p — 1) if a > p + 1 ~ w, unless a > p — r — s. Hence we get: 

p+i-ij 

H{r,s,t;p-l) = - J2 



^p-r\ Ba (p+l-r-s-a\ Sp+i_^_a 



a J p — r \ t Jp+1 — r — s — a 

p-r\ Ba (2p - r - s - a\ B2p-w-a 



a J p ~ T \ t J2p — r — s~a 

r + a\ Ba fw + a-l\ i3p+i-n,-a 



E 

a— max{p+l — r — s,p+2 — It;} 
p-\-l — w 

^ ' \ a ) r + a\ t J w + a-1 
p-i-r (30) 

/ ^^a+t {''" + "■\ _Bafw + a\ B2p — w — a 

, ^ } \a)r + a\t)w + a' 

a— inax{p+l— r — s,p+2— uij 

When w is odd we recover Thm. 13.51 when p > w + 1 by noticing that there are only two nontrivial 
terms in H3U|) corresponding to a = 1 and a = p — w. 

3.4 Some remarkable cases of length 3 with even weight 

Applying (|24|1 to (si, S2) ~ (1, 3) together with shuffle product H{1, l;p — 1) • H{2;p — 1) we get 

i?(2,l,l;p-l) + H(l,2,l;p-l) + i/(l,l,2;p-l) = ^Sp_5 (mod/). (31) 

By Lemma I3. 31 we can even see that 

i?(2,l,l;p-l) = F(l,l,2;p-l) (modp). (32) 

However, is it true that in fact all of these sums are congruent to mod p7 Now from (|3U|I we can 
compute easily that modulo p 

p-3 

i/(l,2,l;p-l)=^B,Sp_3-a, (33) 

a=0 

H{l,l,2;p^l)=-J2^BaBp^3-a, (34) 
I + a 



i7(2,l,l;p-l)=^^^B„i?p_3-a. (35) 

Set 



2 

a=0 



p— 3 

A := BaBp-3-a, B :— ^ aBaBp^^-a- 
Then from (EJ, (EH to ^23) 

A/2 = (modp). 

Moreover, from ^ and ^ 

-{A + B/2) = A/2 + B/2 (modp). 

Consequently we have 

Corollary 3.7. For every prime p > 7 we get 

p— 3 P^3 

BaBp^s.-a = ^ aBaBp^z-a = (mod p). 

a=0 a=0 

Therefore we have 

i/(l,2,l;p-l) = i/(l,l,2;p-l) = H(2,l,l;p-l) = (mod p). (36) 
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Using Maple we further find the following very stimulating example 

2234416196881673576349577192603 _ 2n 
(1,2,1; 36) - 1151149136943530805554073600000 = " (mod 37 ). 

Proposition 3.8. For all prime p > 7 we have 

H{l,2,l;p-l)=-^pBp-5 (modp2), (37) 

i/(2,l,l;p-l)= ^pB.p-5 (modp2), (38) 

H{l,l,2:p-l)= (modp2). (39) 

And 

S{l,2,l;p-l)=-^pBp^5 (modp2), (40) 

5(2,l,l;p-l) = ypBp_5 (modp2), (41) 

^(l,l,2;p-l)= (modp2). (42) 

Proof. Omitting H{- ■ ■ ■,p - I) wc let A = H{1, 2, 1), B = H{2, 1, 1) and C = H{1, 2, 2). Equation 
says that 

A + B + C=^pBp-5 (modp^). (43) 



Now by the shuffie relations 

H{l)H{2, l) = A + 2B + H{3, 1) + H{2, 2), (44) 

From Thm.rOwe get 



3 

lb 

Hence it suffice to show (|37|) . From (|12|l we see that 



A + 2B = —pBp_5 (mod p^) (45) 



S'(l, 2, 1) = H{1, 2, 1) + (3, 1) + H{1, 3) + H{A) 

= i?(l,2,l) + i?(l)-iJ(3) = iJ(l,2,l) (modp^) (46) 
S{2,2)=H{2,2) + H{A) = -H{2,2) + H{2f = -H{2,2) (mod p^). (47) 



Now because (1, 2, 1)* = (2, 2) by Thm. 12111 we have 

-S'(2,2) = S'(l,2,l)+pfi/(l,2,l,l) + i/(3,l,l) + iJ(l,3,l) + i?(4,l)) (modp^). (48) 



By Cor. l3.4l iJ(l. 3, 1) = (mod p). So using expressions H46|) and (|47|l we can simplify the preceding 
congruence to 

iJ(2,2) = F(l,2,l)+p(ff(l,2,l,l) + F(3,l,l) + F(4,l)) (mod p^). (49) 

Now that (1, 2, 1, 1)* = (2, 3) we find from congruence (|10|l that modulo p 

~ S{2, 3) = S{1, 2, 1, 1) = H{1, 2, 1, 1) + H{3, 1, 1) + H{1, 3, 1) 

+ H{1, 2, 2) + H{4, 1) + H{3, 2) + H(l, 4) + H{5) (mod p). 

Note that 

i7(l, 3, 1) = 0, H{A, 1) + H{1, 4) = i/(l) • H{A) - H{5) = (mod p). 
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So ((T^ implies that 

-i7(2,3) = -S'(2,3) = iJ(l,2,l,l) + ff(3,l,l) + F(l,2,2) + ff(3,2) (modp). 

Namely 

2, 1,1) = - i/(2, 3) - His, 2) - i/(3, 1, 1) - Hil, 2, 2) (mod p) 
= H{5) ~ H{3,1,1) - H{1,2,2) (mod p) 

= -iJ(3,l,l)-i7(l,2,2) (modp). 

Plugging this into (|49|l we see that 

H{2, 2) = H{1, 2, 1) + p(i?(4, 1) - iJ(l, 2, 2)) (mod 
Using Thm. E!T^ and Thm. I!01 we can now compute easily that 

i/(2,2) = -^pBp_5 (mod/), i7(4,l) = -Bp_5, F(l, 2, 2) = -^Bp_5 (mod p). 

These lead to congruence H37f) . Then congruence H4Uf) follows from 1)46(1 . 

We now can solve H43|) and (|45|l to get congruences H38() and H39|l . Finally, (|13|l yields 

5(2, 1, 1) = H{1, 1, 2) - (l)i/(l, 2) - H{1, 1)H{2) + H{lfH{2) = ^pBp-5 (mod p^), 
5(1, 1, 2) = Hi2, 1, 1) - Hi2)Hil, 1) - iJ(2, l)i7(l) + H{2)H{lf = ^-pB^^^ (mod p^). 

We have completed the proof of the proposition. □ 

By going through the proof of Thm. rTHl or the proof leading to 130(1 we can obtain the following 
result. We leave the details of the proof of it to the interested readers. 

Proposition 3.9. For all prime p > 17 we have 

77(4, 3, 5; p- 1) = ff(5,3,4;p- 1) = (modp). 

In [38( . by using the shuffle relations and Hoffman's convolution we will study the mod p struc- 
ture of the multiple harmonic sums for lower weights. In particular, we will prove Prop. 13.91 and 
congruences like 

5(2, 3, 2, 3, 2;p - 1) = 5(2, 3, 3, 2, 2;p - 1) = (mod p). 

Problem 3.10. Numerical evidence shows that if (r, s, t) ^ (5, 3, 4), (4, 3, 5) and r + s + t>Q\s even 
then the density of primes p such that p \ H{r, s,t]p — \) among all primes is always 1; however, 
there is always p such that p\H{r, s,t;p — 1). Can one generalize the formula in Thm. 13.51 to prove 
this? 

3.5 Some congruences modulo prime squares 

We know from Thm. 13.51 that H{r, s,r;p — 1) = (mod p) if s is an odd number and p > 2r + s. 
Modulo p^ we have (using substitution of indices k —> p ~ k) 

H{r, s,r;p — 1) = —H{r, s,r;p — 1) — p^2rH{r + 1, s,r;p — 1) + sH{r, s + 1, r;p — 1)] . 

Hence 

iJ(r, s, r;p — 1) = — p[riJ(r + 1, s, r;p — 1) + -i7(r, s + 1, r;p — 1)] (mod p^). (50) 
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Proposition 3.11. For all prime p > 5 we have 

i7(l,3, 1) =5(1,3, 1) = (mod p^) (51) 

77(2, l,2;p-l)= 5(2, l,2;p-l)^-ii32_3p (modp^). (52) 
For all prime p > 7 set bg,{p) — (55(6, 1, l;p — 1) + Bp^^^Bp^^) /2 then we have 

i?(l,5,l;p-l) =5(l,5,l;p-l) = 68(p)p (modp^) (53) 
i?(2,3,2;p-l)=5(2,3,2;p-l)= 468(p)p (mod/), (54) 
ff(3,l,3;p-l)=5(3,l,3;p-l) = 68(p)p (modp^). (55) 

Proof. Notice that 

5(r, s, r;p — 1) = H{r, s, r;p — 1) + H{r + s)H{r) = -ff (r, s,r;p ~ 1) (mod p). 
From H50() we have 



2i/(l,3,l;p- 


1)^ 


-p[2iJ(2,3,l;p- 


1)- 


h3i/(l,4,l;p 


-D] 


(mod 


2i/(2,l,2;p- 


1)^ 


-p[4i/(3,l,2;p- 


1)- 


Fff(2,2,2;p- 


D] 


(mod p^), 


2if(l,5,l;p- 


1)^ 


-p[2i?(2,5,l;p- 


1)- 


h5i/(l,6,l;p 


-D] 


(mod p^ 


2F(2,3,2;p- 


1)^ 


-p[4i?(3,3,2;p- 


1)- 


h3i?(2,4,2;p 


-D] 


(mod 


2i?(3,l,3;p- 


1)^ 


-p[6i/(4,l,3;p- 


1)- 


hF(3,2,3;p- 


D] 


(mod p^). 



By gni Thm. 7.2] we know that 

5(2,3, l;p- 1) =35(4, l,l;p- 1) (modp), 
5(1,4, l;p- 1) = - 25(4, l,l;p- 1) (modp), 

5(3,l,2;p-l)= -5(4,l,l;p-l) = iB2_3 (mod p) 

which yield 1)51(1 and ((52|l . Similarly, ((53(1 follows from 26, Thm. 7.4] because 

25(2,5, l;p- 1) = -25(5,2,l;p- 1) - 25(5, 1, 2;p - 1) (modp), 

= 55(6, 1,1; p-l)-Bp_5Sp_3 (modp), 
5(l,6,l;p- 1) = -25(6,l,l;p- 1) (modp). 

Also from 26, Thm. 7.4] 

5(3,3,2;p-l) =4fe8(p) (modp), 

5(2, 4, 2;p - 1) = - 25(4, 2, 2;p - 1) = -25(3, 3, 2;p - 1) (mod p), 

5(3, 2, 3;p - 1) = - 25(3, 3, 2;p - 1) (mod p), 

5(4,l,3;p-l) =fo8(p) (modp). 
These lead to the last two congruence of the proposition immmediately. □ 

By similar argument we can compute H{r, s, r) (mod p^) for odd s if we know the values 5(s;p— 1) 
(mod p) for all s of length 3 and weight 2r + s + 1. When we apply this argument to s = (1, 3, 1, 3) 
we find 

Proposition 3.12. The following two congruences are equivalent: 
(i) For all primes p > 8 

31 

7J(l,3,l,3;p- 1) = -— pBp_9 (mod/). 



(ii) For all primes p > 8 

1^3 1889 
54 648 

We have verified the congruences in the proposition for all primes p such that 10 < p < 2000. 



5(6, 1, 1, l;p - 1) ^ - -^Bp^9 (mod p). 
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3.6 Some congruences of Bernoulli numbers 
From Cor. 13.71 we see that for every prime p > 7 we have 

p— 3 

^ BaBp_3_a = ^ aBaBp^s^a = (mod p). 

a=0 a=0 

Can we generahze this? The answer turns out to be affirmative. 
Proposition 3.13. For every prime p > 9 we have 

p-5 2 

^BaBp_5_a = --S2_3 (mod p) . 

Proof. By (I30|) we have for any even number n 

p— n— 1 p— 2 

H{l,n,l;p-l)= J2 (-l)"S„Bp_„_i_„ + ^ (-l)"B,B2p-„-2^a (56) 

a— a—p—n 

Taking n — A and comparing with |26l Thm. 7.2] we get 

p-5 ^ 

^i3„Sp_5-a + = ^(1,4, 1) = 5(1,4, 1) = -i?2_3 (mod p). 

a=0 

This proves the proposition. □ 

The foUowing result is straight-forward. 
Proposition 3.14. For all positive number n and prime p > n + 3 we have 

i/(l,l,l,n;p-l) = -(-l)"i/(7i,l,l,l;p-l) = (-ir5(l,l,l,77;p-l) = -5(n,l,l,l;p-l) 

p — 2p — 71 — 1 , 7\ / ; \ j-> 7-> 7-> 

'a + b\ /a + b + n\ BaB^Bp^n-i-a-b 



a=0 fc=0 



^ ' \ h l\ r, I {a + l){a + b+l) ^ ^' 



3.7 Multiple harmonic sums of arbitrary length 

To prove the main resuh in this section let us recall the Bernoulli polynomial {x) which is defined 
by the following generating function 



^ j-m 

= y B^{x)—. 

1 m! 



te 

m=0 



These polynomials satisfy (see [271 p. 248]): 



kx^~^ (57) 

B'„^{x) ^mBrn-i{x) Vm > 1. (58) 
Lemma 3.15. For any even integer n > and prime p > 3n + 7 we have 

a,b>0, a+b=p-3n-3 ^ ^ ^ ' 

y A('Tt)^('T:)-0 (modp). (60) 

^-^ p ~ b \n + 1 / p — a \n + 1 / 

a,6>0, a+6=p-3ji-3 ^ ^ ^ ' 
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Proof. Throughout the proof all congruences are modulo p. First we set 

_Bh_fp-b\^a_fp- a 
^ p - b\n + 1 J p - a\n + 2 

V _Bb_/'p-b\^a_/'p- a 

^-^ n — 6 In + 1 / p — a \n + 1 



We have 



n + 2 -p — h\n^ \) p — a\n+ \ 

a.b>0, a+b—p—3n~3 



Now exchange the index a and 6 in A we get 

Bb f p~b\ Ba (p~ a 



p — b\n -\- 2 I p — a\n \ 

271 + 2 + Q Bb fp-b\ Bg fp-a\ 
^ n + 2 p-b\n+\j p- a\.n+lj' ' 

a,b>0, a+6=p-3n-3 ^ ^ ^ ^ ^ 



E ^ BbfP~b\^h^a\^n^ 

^ 71-1^9 Tl— /)\n-4-l/T)— 77-1-9 ^ ' 



Adding (jHH) to ^ we find 

^ 1,1,- 

n + 2 p — b\n + Ij p — a\n + IJ n + 2 

a,6>0, a+b— p~3n — 3 

Now for any nonnegative integer m and h such that 6 + m + 1 < p we have 
We get 

a,b>0, a+b=p-3n-3 ^ ^ ^ ^ 

because all the terms with odd index a are zero. Let us denote by Coeff[i, f{x)] the coefficient of 
in the polynomial f{x). Then (|65|l means that 



B = n^Coeff[p + n + 1, B;_„_i(2;)] 

9 

^ Cocff[p + n, (i?^_„_i(a:)) ] 

p + n + 1 ^ 



n + 

by (jsHJ. From ifFfI) this yields 



YCoeff[p + n, 2(p - n - f )Bp_„_i(a;)-Bp_„_2(x)] 



p — n — 1\ fp — n — 2 



a,b>0, a+b— 377 — 3 

_ 2n Sfc f p — b\ Ba fp — a 



n + 1 ^-^ p — b\n + 1 / p — a\n + 2 

a,b>0,a+b=p-3n-3 ^ ^ ^ ^ ^ 

by Therefore 

B = (66) 
n + f ^ ^ 

Now the lemma follows easily by (|63f) and H66|) . □ 
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Theorem 3.16. Let p be a prime and s S N'. Assume p > |s| + 2. Then 

H{s;p~l) = S{s;p-l)=0 (mod p) 

if s has the following forms: 

1. s ^ (l™, 2, 1") for TO, 71 > and m + n is even. 

2. s= (l", 2, 2, w/iere n > 2 is even. 

3. s = (l"+\2,l"-i,2,l") u;/ieren>0. 

^. s = (l", 2, 1", 2, 1") where n>2 is even. 

Furthermore, in the first and last cases we also have S(s;p — 1) = (mod p). 

Proof. We omit p — 1 in S'(— — 1) and H(—]p — 1) and assume all congruences are modulo p 
throughout this proof. 

Let s — (1™,2,1"). For future reference we first allow m and n to have different parity. 
Clearly we have s* = (to + 1, n + 1). So |s| = |s*| = to + n + 2, Z(s) = to + n + 1, and ^(s*) — 2. By 
123 Thm. 6.7] we have 

S'(s*) = -S'(s). (67) 

Since |s*| - ^(s*) = to + n we have by 

(-l)"+"5(s*) = (-l)l"'l^'("*)^(s*) = ~H{s*) + H{m + 1) • H{n + 1) = -H{s*) (68) 
by Thm. If .31 We now apply (|f 3|1 to s and get 

-5(s)^ ^ i-lYllHis,). (69) 
Now if Z > 2 then one of Sj = for some positive d so that H{sj) = 0. Thus 

^(S) EE H{S). 

Combining this with (|67|l and H68|) we get 

H{S) EE (-f)™+"i/(s*). (70) 

When to + n is even we get 

S{S) EE H{S) EE H{S*) EE 0. 

by Thm. 13.21 which proves case . 

131 Let s = (l",2,l"-\2,f"+i) where n > 2 is even. Then s* = {n + l,n + l,n + 2). So 
|s| = |s*| = 3n + 4, ^(s) = 3n + 2, and ^(s*) = 3. Since n is even we have by applying (|f 3|l to s* 

S{s*) EE H{s*) - H{2n + 2,n + 2)- H{n + f , 2n + 3) + H{3n + 4) = H{s*) (71) 

by Thm. Inland Thm. [TTSl Applying (|13|) to s and using the fact that H{1'^) = for any d we have 

ri-l 

S{s) EE (-l)l"l-'(")S'(s) = - H{s) + J2 H{V\ 2, 1'') • H{V'-^-'', 2, 1"+^) 

n-l 

EE - H{s) + ^H{n+l,a+l) ■ H{n -a,n + 2) 

by GOJ). Hence by (E3 and GH) 

n-l 

H{s) EE H{n + l,n+l,n + 2) + H{n + 1, a + 1) • H{n -a,n + 2). (72) 

a=0 
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We know that for all j,k < p we have 

p-2 
a=0 



(mod p) ii j k, 
— 1 (mod p) if j = k. 



It follows that 

p-2 p-2 ^ ^ 

EH(n + l,a+l) ■ H(n- a,n + 2) V V . . -, r-pr 

a^O a^O l<z<j<p l</c</<p 



E 



l<i<j=A:<i<p 

= - i?(n + l,n + 1,71 + 2). 

Together with (|72(l we see that 

p-2 

H{s) = -Y^ H{n + 1, a + 1) ■ H{n - a, n + 2) 

a—n 

{-lYBp-n-2-a fn+2 + a\ Ba^2n-1 + 2n + 1 - a 

~ ^ n + 2 + a V «+l /p + 2n+l-aV n + 2 

a— 2n+l \ / -I \ 

by Thm. 13.11 Under substitution a ^ 2n + 1 + a we get: 

^(s)^'EVir("":!|")^^^(^T:)^-0 (mod,) (73) 

a— ^ ' \ / ^ 

by H59|l . This combined with H71I) completes the proof of case (|2J). 
Q It follows from (0) by taking s. 

I0J When 71 = or n is odd this follows from Thm. [TT^ and Lemma [3.31 respectively. When 
71 > 2 is even the proof is almost the same as that of case ^ except at the end one need resort to 
l|Hn|. The congruence for S follows from the fact that (l", 2, 1", 2, 1")* = (71 + 1, ti + 2, ii + 1) and 
therefore 

S'(s*) = i/(s*) = 

by Cor. EH □ 

Remark 3.17. Note that in cases Q and we usually have iS'(s) ^ (mod p). For example, in 
case we have S'(s) = -S'(s*) = -i?(s*) = H{n + 1, n + 1, n + 2) (mod p) by (jZHl. We know that 
iJ(3, 3, 4; 12) = 8 (mod 13), i7(5, 5, 6; 18) = 15 (mod 19) and H[l , 7, 8; 28) = 26 (mod 29). 

Theorem 3.18. Let s = {r, s}" /or so77ie r, s > 1 and p > \s\ be a prime. Then 

H{s;p-l) = S{s;p-l) = (mod p) 
if either (i) n — 1,2, both r and s are even, or (ii) n is any positive integer, both r and s are odd. 

Remark 3.19. When 71 1 this is Thm. EH When n = 2, this has been confirmed by Hoffman 
(see the remarks after Thm. 6.3]). 

Proof. By the above remark we may assume r and s are odd and proceed by induction on 71. In the 
following we will drop p — I again. By the shuffle relations and equation (|12|l it is straightforward to 
verify that 

S{{r,sr)^ J2 Hit) 

t^{r,s}" 

= H{{r,sr) + H{r + s)- ^ H{t) 

t^{r,s}"-i 

= 77({r,s}") + H(r + s) • S'({r,s}"-i) = H({r,s}") (mod p) (74) 
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by induction assumption. On the other hand, from we find that 

5({r,sn^ J2 (-l)'ff(si)---i?(sO (modp). 

Now if / > f then Si — {r,s}'^ for some d < n in which case H{si) = (mod p) by induction 
assumption, or else, Si = {r, s}'' U {r}, in which case Si is a pahndrome of odd weight and hence 
ff(si) = (mod p) by Cor. 13.41 Consequently 

S{{r,.sr)^-H{{r,sr) (modp). 

Together with H74|l it shows that H{{r, s}") = S{s;p—1) = (mod p) and the theorem is proved. □ 

When both r and s are even but n > 2 the theorem does not hold in general. For example, 

i/({2,4}3;22) = 21 (mod 23), i7({2, 4}^; 28) = 20 (mod 29). 

3.8 Some conjectures in general cases 

When the length I > 4 numerical evidence up to length 40 shows the following conjecture is true. 
Conjecture 3.20. Let s G and p > |s| be a prime. Then 

H{s;p—1)=0 (mod p) 

if s has one of the following forms: 

1. s ~ (1 1™, 2, 1", 2}^, 1™, 2, 1") for q^m,n > 0, where either (i) q is odd, or (ii) q is even and 
m + n is even. 

2. ({2}", {3, {2}"}") form,n>0. 

3. (l,{2}'", {1,{2}™+1}",1,{2}'",1) form,n>0 and n is even. 

We conclude our paper by 

Problem 3.21. Are there any other s satisfying Wolstenholme's type theorem besides those we 
listed in the paper? More generally is it possible to find a formula similar to Thm. 13.51 for arbitrary 
s? 

Appendix: Distribution of Irregular Primes 

Tabled in the paper and the following Table [21 give us some evidence to Conjecture 12. 171 

In Table 121 we count the first 30,000 irregular primes with irregular index ip = 1, first 15,000 
irregular primes with index 2 (producing 30,000 irregular pairs), and all the irregular primes < 
12,000,000 with index 3 (producing 3 x 9824 = 29472 irregular pairs). We denote by N{k,m) the 
number of irregular pairs (p, t) satisfying p — t = k (mod 3) in the top m irregular pairs, < A: < 2, 
and by P{k,m)the percentage of such pairs. We put a subscript i for irregular primes of index i. 

Between 3 and 12 million there are only 1282 irregular primes with irregular index 4 (producing 
5128 irregular pairs), and 127 irregular primes with irregular index 5 (producing 635 irregular pairs). 
We provide the data for them in Table 13 and Table 01 There are 13 irregular primes < 12,000,000 
with irregular index 6, producing 78 pairs. For them we have ^^6(0, 78) — 25, Ne{l,78) = 27, 
and Ne(2, 78) = 26. There are merely 4 irregular primes in the same range with irregular index 7, 
producing 28 pairs. For these pairs: iV7(0,28) = 10, Nril,28) = 9, and ^7(2,28) = 9, which is the 
best we can hope. 
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ip — 1 , m 


3,000 


6,000 


9,000 


12,000 


15,000 


18,000 


21,000 


24,000 


27,000 


30,000 


iVi(0, m) 


979 


1968 


2954 


4018 


5001 


5993 


6972 


7973 


8968 


9993 






oz.oL) 


OO OO 

OZ.OZ 


O O AO 


O O O /I 


O O OA 

oo.zy 


O O OA 

00. zL) 


O O OO 

OO.ZZ 


O O OO 

OO.ZZ 


o o o 
OO.O 


AT /I 

iV 1 ( 1 , m j 


lUiD 


ZUZD 


o A /I n 


/I A /I O 


c Aon 
51)09 


o(J ib 


TAA C 


Olio 

olio 


AAAA 
91)91) 


1 A A C C 


_ri (1, mj 


33.87 


33.77 


o o o 

33.88 


O O o 

33.68 


*) o e A 

33.59 


O O 

33.75 


O O i^A 

33.79 


o o o o 

33.82 


33.67 


33.52 


AT /O 

iVi(z, mj 


iUUo 


2006 


299 ( 


OA /I A 

3940 


A A^ A 

4960 


C AOO 

5932 


AO O 

6933 


^A AA 

/yo9 


O A /i O 

8942 


A A C O 

9952 


Pi(2,m) 


OO. oU 


33.43 


O O OA 

33.30 


OO o o 

32.83 


OO A7 

33.0/ 


OO A^; 
32.96 


O O A 1 

33.01 


OO AC 

32.95 


O O 1 o 

33.12 


O O 1 

33.1 / 


Zp = 2, TO 


3,000 


6,000 


9,000 


12,000 


15,000 


18,000 


21,000 


24,000 


27,000 


30,000 


A^2(0, to) 


981 


2029 


3033 


4049 


5013 


6024 


7064 


8062 


9085 


10096 


Jr 2 [yj- III) 


,>2. 1 U 


33.62 




66. 1 4 


O O 1 o 

00.42 


O O 1 

33.4 < 


00.04 


33. -I)!) 


OO 1' r 


o o 1' r 


iV2(i, mj 


993 


1977 


3033 


4005 


5001 


6012 


AO 1 

6981 


7962 


8927 


9910 


P2{l,m) 


33.10 


32.95 


O O 'TCX 

33.70 


o o o o 

33.38 


O O OA 

33.34 


O O A C\ 

33.40 


O O O /t 

33.24 


O O 1 o 

33.18 


O O A/^ 

33.06 


O O AO 

33.03 


7V2(2, mj 


1026 


1 AA /I 

1994 


ono A 

2934 


3946 


4986 


c: A/^ A 

5964 


6955 


TAT/^ 

Tdlb 


O AO O 

8988 


AAA A 

9994 


ii2(2, TOj 


O A OA 

34.20 


33.23 


OO 

32.00 


OO o o 

32.88 


O O O /i 

33.24 


O O 1 o 

33.13 


O O 1 o 

33.12 


OO OO 

66.16 


O O OA 

33.29 


O O O 1 

33.31 


ip = ?>,m 


3,000 


6,000 


9,000 


12,000 


15,000 


18,000 


21,000 


24,000 


27,000 


29,472 




997 


2027 


3065 


4102 


5124 


6091 


7087 


8067 


9046 


9889 


P3(0,to) 


33.23 


33.78 


34.05 


34.18 


34.16 


33.84 


33.74 


33.61 


33.50 


33.55 


iV3(l,m) 


1001 


1994 


2968 


3955 


4944 


5956 


6921 


7943 


8974 


9812 


P3(1,to) 


33.36 


33.23 


32.97 


32.96 


32.96 


33.09 


32.95 


33.09 


33.23 


33.29 


A^3(2,m) 


1002 


1979 


2967 


3943 


4932 


5953 


6992 


7990 


8980 


9771 


P3(2,m) 


33.40 


32.98 


32.96 


32.86 


32.88 


33.07 


33.29 


33.29 


33.26 


33.15 



Table 2: Distribution of p — t (mod 3) for irregular pairs {p, t) with ip = 1, 2, 3. 



ip = 4, TO 


600 


1,200 


1,800 


2,400 


3,000 


3,600 


4,200 


4,800 


5,128 


7V4(0,to) 


196 


388 


569 


769 


948 


1161 


1341 


1537 


1645 


P4(0,to) 


32.67 


32.33 


31.61 


32.04 


31.60 


32.25 


31.93 


32.02 


32.05 


A^4(1,to) 


202 


390 


615 


818 


1032 


1217 


1428 


1624 


1734 


P4(1,to) 


33.67 


32.50 


34.17 


34.08 


34.40 


33.81 


34.00 


33.83 


33.79 


^4(2, to) 


202 


422 


616 


813 


1020 


1222 


1431 


1639 


1749 


P4(2,m) 


33.67 


35.17 


34.22 


33.88 


34.00 


33.94 


34.07 


34.15 


34.08 



Table 3: Distribution of p — t (mod 3) for irregular pairs {p, t) with ip = 4. 
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ip = 5,m 


100 


200 


300 


400 


500 


600 


635 


N'^i O, m) 


31 


71 


111 


145 


179 


218 


230 


^5(0. //)) 


31.00 


35.")() 


37.00 


30. 2.") 


35.80 


3G.33 


35.94 


iV5(l,m) 


33 


58 


87 


115 


154 


190 


200 


P5(l,m) 


33.00 


29.00 


29.00 


28.75 


30.80 


31.67 


31.25 




36 


71 


102 


140 


167 


192 


205 


P5(2,m) 


36.00 


35.50 


34.00 


35.00 


33.40 


32.00 


32.03 



Table 4: Distribution of p — t (mod 3) for irregular pairs {p, t) with ip = 5. 
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